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ABSTRACT 

We report on the identification of a third, new ultraluminous X-ray source (ULX) in the starburst 
galaxy M82. Previously, the source was observed at fluxes consistent with the high state of Galactic 
black hole binaries (BHBs). We observe fluxes up to (6.5 ± 0.3) x 10^^ ergs s~^ in the ultraluminous 
regime. When the source is not in the the low/hard state, spectral fitting using a multicolor disk 
model shows that the disk luminosity varies as the disk inner temperature raised to the power 4.8±0.9, 
consistent with the behavior of Galactic BHBs in the thermal dominant state. Fitting the spectrum 
with a multicolor disk model with general relativistic corrections suggests that the source harbors a 
rapidly spinning black hole with a mass less than 100 solar masses. A soft excess was found in the 
source spectrum that could be blackbody emission from a photosphere created by a massive outflow. 
The source also showed soft dips during a flare. 

Subject headings: accretion, accretion disks — black hole physics — X-rays: individual (CXOM82 
J095546.6-F694041, X37.8-h54) 



1. INTRODUCTION 

Ultraluminous X-ray sources (ULXs) are nonnuclear, 
point-like X-ray sources in external galaxies with an ap- 
parent luminosity over the Eddington limit of a stellar- 
mass black hole, which is typically 3 x 10'^^ ergs s~^ for 
a 20 Mq black hole. Those with variability on time 
scales from minutes to months are accreting compact 
objects. The nature of ULXs is still unclear. If the 
emission is (roughly) isotropic and under the Edding- 
ton limit, they may harbor black holes up to 10'^ Mq, 
belonging to the missing population of interme diate 
mass black holes (IMBHs; 'Colbert & Mushotzkv' IW 
Makishima et al. 2000; Kaaret et al. 2001; Farrcll et al 
20091). However, they could also be stellar- mass 



black holes, with b eamed or supe r-Eddi ngton emission 
(iKing et all 120011: iKording et al.l 120021: iWatarai et all 
l200lHBegelmanll2002t ). 

ULXs with relatively low luminosities (^ 3 — 6 x 
Iq39 gj.gg g-i-j j^g^y |-jg jjiorc Hkcly to be stellar-mass 
black holes of about 10 Mq like Galac tic BHBs or up 
to a c ouple tens of Mq like IC 10 X -1 (jPrestwich et al.l 
[20071 : [Silverman fc Filippenkd [2008[ ). Galactic BHBs 
have been observed in a few cases with luminositie s up 
to 6 X 10^^ ergs s'^ (|McClintock fc Reniillardll2006[) . al- 
though there ar e large uncertainties on their distances. 
[Gladstone et al.l (|2009f l proposed that most ULXs are 
stellar-mass black holes in an 'ultraluminous state' with 
super-Eddington accretion based on analysis of high 
quality X-ray spectra. A major difference between the 
two populations is that Galactic BHBs spend most of 
their lifetime in the quiescent or low state, while most 
ULXs have persistent emission and have been active since 
their discovery at time scales up to decades in some cases. 
The non-transient behavior of the majority of ULXs im- 
plies that they do not contain blac k holes s ignificantly 
more massive than 10^ Mq (jKalogera et al.|[2004i) . 
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spite these hypotheses, the connection between ULXs 
and Galactic BHBs has not been well addressed. If they 
contain black holes of similar mass, we may expect to see 
transitions from the normal states to the ultraluminous 
state in an individual source. 

Multicolor blackbody emission from an optically thick 
accretion disk could shed light on the nature of the cen- 
tral black hole. At a fixed fraction of Eddington luminos- 
ity, the disk inner temperature scales with th e black hole 
mass to the —1/4 power (fMakishim a et all [2000'). indi- 
cating that accretion disks around IMBHs are cooler and 
brighter than those around stellar-mass black holes. Soft 
excesses are detected in ULX s and can be modele d as cool 
disk emission from IMBHs (jKaaret et al.l [20031) . How- 
ever, the cool disk model is not the unique interpretation 
of the soft excess (for a brief review see the introduction 
of Feng & Kaaret 2009). Massive outflows due to near or 
super-Eddington accretion may pr oduce blackbody emis- 
sion at the similar temperature (IKing fc Poun ds l2003|; 
[Km3[200l iBegelman et "aI|[20Q6l: [Poutanen et al. 2007). 

The disk model can be tested by determining if the 
disk luminosity varies with the 4th power of the in- 
ner disk temperature, i.e. idisk oc T^. This relation 
has been robustly verified for stellar-mass black holes 
when the disk emission dominates the X-ray spectrum 
([Gierhhski fc Dond 120041 ). Although the X-ray spectra 
of so me ULXs can be modeled by dominant d isk emis- 
sion (jStobbart et al.l 120061 : iWinter et ahl 12006") . the ob- 
served luminosity variation and number of observations 
have been inadequate for proper testing of the Ldisk oc T^ 
relation. A possible L ex T^ relation was found in NGC 
5204 X-1, but the disk emission was not the dominant 
component in the spectrum a nd a la rge correction for 
Comptonization was required (jFeng fc Kaaret 2009) . 

M82 is a nearby starburst galaxy in wh i ch two ULXs 
have been found (iMatsumoto et al.l 12 0011: iK aaret et al.l 
[200lt iFeng fc Kaareti I2007bl : IKonTet al. 200^. It has 
been extensively studied in X-rays with the Chandra 
X-ray observatory . In our recent Chandra observa- 
tion on 2008 October 4, the X-ray source CXOM82 
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Table 1 

Chandra ACIS Observations of M82 used for spectral analysis in 
the paper. 



ObsID 



361 
1302 

379 

2933 

6097 

5644 

6361 

8190 

10027 

10025 

10026 



Date 


Instrument 


1999-09-20 


ACIS-I 


1999-09-20 


ACIS-I 


2000-03-11 


ACIS-I 


2002-06-18 


ACIS-S 


2005-02-04 


ACIS-S 


2005-08-17 


ACIS-S 


2005-08-18 


ACIS-S 


2007-06-02 


ACIS-S 


2008-10-04 


ACIS-S 


2009-04-17 


ACIS-S 


2009-04-29 


ACIS-S 



Exposure 



33.3 
15.5 
8.9 
18.0 
52.8 
68.1 
17.5 
52.8 
18.3 
17.4 
16.9 



J095546.6-f694041 was 5-10 times brighter than its nor- 
mal flux level and was identified as the third ULX in 
M82. Accord i ng to the naming convention suggested by 
iKaaret et alJ (|2001l ). we refer to this source as X37. 8-1-54. 
In this paper, we investigated its spectral and timing 
behavior using all available Chandra observations, and 
discussed its possible nature. The distance to the h ost 
galaxy is adopted as 3.63 Mpc (jFreedman et al.iri994[ ). 

2. OBSERVATIONS AND DATA ANALYSIS 

So far there have been fourteen observations of M82 
with the Chandra Advanced CCD Imaging Spectrome- 
ter (ACIS), including the most recent three that we pro- 
posed to monitor the brightest ULX in it (the results 
will be reported separately). Among them, eleven obser- 
vations with an exposure longer than 5 ks are used for 
spectral analysis (Table [1]). The other three short ACIS 
observations (Obs ID 378, 380-1, and 380-2; 380 has two 
observations) as well as the four High Resolution Cam- 
era (HRC) observations (Obs ID 1411-1, 1411-2, 8505 
and 8189; 1411 has two observations) are examined for 
long term variability. All data are processed using CIAO 
4.1.2 with CALDB 4.1.3; new bad pixel files and level 2 
events files are recreated when necessary. 

Two images of the central region of M82 in the 0.3- 
8 keV energy range taken from Chandra ACIS are dis- 
played in Figure [T] The top panel shows the image from 
observation 10027 when X37. 8-1-54 was the brightest ever 
seen, and the bottom one is from observation 5644 in 
which the second brightest fiux was recorded. The source 
in 10027 spreads on multiple pixels because it lies at 
a large off-axis angle. Aligning the two images using 
isolated point-like sources suggests that the brightening 
in observation 10027 around the source region is indeed 
from X37.8-f 54. A faint object to the south of X37.8-t-54 
is seen in observation 5644, and is partly blended with 
X37.8-f 54 only in observation 10027. It displays a vari- 
ability by a factor of 2. For observation 10027, we sub- 
tracted the contribution of this dim source assuming a 
spectrum from its highest state and found no significant 
change in the spectral fitting for X37. 8-1-54, indicative of 
negligible contamination. 

The energy spectra of X37. 8-1-54 are extracted from 
events in a 3 tr elliptical region around the source found 
by wavdetect. The background is subtracted from a 
nearby circular region with a radius of about 4 arcsec 
on the same CCD and off the readout column of the 
source. Background lightcurves are checked using 10-15 
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Figure 1. Chandra ACIS image of the central part of M82 in 
the energy range of 0.3-8 keV from observations 10027 (top) and 
5644 (bottom), respectively, where X37. 8-1-54 was in outbursts. The 
pluses indicate the three ULXs in M82. 



keV events on the whole CCD; intervals with background 
flares were excluded in observations 5644 and 6361. Spec- 
tral fitting was done using XSPEC 12.5.0 and the errors 
are quoted at 68% confidence level if not specified other- 
wise. Observations 361 and 1302 were performed in the 
same day and the spectral shapes are consistent. These 
two spectra are fitted together to improve the signal 
to noise ratio. The Galactic absorption column density 
along the line of sigh t to the source is 5.05 x 10^° cm~^ 
(Kalb erla et al.|[2005D . and is set as the lower boundary of 
the total absorption column density during the fits. Us- 
ing the pileupjnap tool available in CIAO, we obtained 
that the count rate on the brightest 3 by 3 pixels is < 0.1 
counts per frame. For the two observations 5644 and 
10027 when the source was brightest, the highest count 
rate in a 3 by 3 pixel island is about 0.03 counts per 
frame due to a large off-axis angle and/or a small read- 
out time. The derived spectral parameters do not change 
within errors when a pileup model is applied. Therefore, 
the effects of pileup are insignificant and, thus, not taken 
into account. 
First, we tried a power-law model subject to interstel- 
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Table 2 

Spectral parameters of X37.8+54 derived from the power-law model subject to interstellar absorption. 



ObsID 


A^H 


r 


JVPL 


/x 


Lx 


xVdof 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


361+1302 


7.71°-; 


2 60+°-" 


9 0+0.3 
^"^-0.3 


0.302t°:°ll 


1 7+0.3 
^•'-0.2 


60.9/54 


379 


3.6ti:^ 


2 0+°-2 
^■^-0.3 


i-8lg:^ 


0.57t°:°^ 


1 i;+0.3 
J^-^-0.2 


6.4/12 


2933 


3 2+°-7 


1 62+°12 


0.92t°:l«3 


n 4^5+0.03 


„ Q7+0.06 


19.9/28 


6097 


6 7+°* 


2 qq+°" 


9 p,+0.2 


9no+°°°8 

0.^00_(,Qgg 


1 fi+°-3 


37.9/50 


5644 


7.0t°J 


2.64«:S^ 


c q+0.3 


„ 701 +0-014 


4 4+0-2 


161.2/149 


6361 


6.4«:? 


2.68tg:}« 


9 q+0.5 


n qq+002 
'^••^^-0.02 


9 9+0.5 


20.3/25 


8190 


4 Q+0-3 


r, 99+0.07 
^■^^-0.07 


1 7c:+0.13 


n4oi+°°i2 


1 or:+0.08 
J^-'^°-0.07 


92.6/82 


10027 


7 c+0.4 
'•^-0.4 


1 S«+006 
^■»»_0.08 


7.4i°:« 


9 99 + 0.05 


0-^-0.3 


122.4/115 


10025 


2.6tS-i 


-, 70+0.14 


0.98t°:ll 


PI 40+0.03 


nq^+0-07 


33.3/26 


10026 


4.4t°;? 


-, 77+0.14 
^■''-0.13 


1.071°;}^ 


p, 4, +0.03 
^•41-0.03 


1 no+°°s 


19.2/23 



Note. 

Col. (4): 



units of 10^^ cm ^. Col. (3): The power- law photon index. 
2 s-l. Col. (5): Observed flux in 0.3-8.0 keV in units of 
Col. (6): Unabsorbed luminosity in 0.3-8.0 keV in units of 10^® ergs s~^ assuming istropic emission at a distance of 



— Col. (1): Observation ID. Col. (2): Absorption column density in 
The power-law normalization at 1 keV in units of 10~* photons cm 



3.63 Mpc. Col. (7): Best-fit x and degrees of freedom. 



Table 3 

Spectral parameters of X37.8+54 fitted with the MCD model or the MCD plus blackbody model subject to interstellar absorption. 



ObsID 

(1) 


TVh 

(2) 


Tin 

(3) 


RinVcOsi Tbb 

(4) (5) 


fx 

(6) 


Ldisk cos i 
(7) 


-f-bb 

(8) 


xVdof 
(9) 


Model: MCD 



379 

2933 

10027 

10025 
10026 



-1.1 

-0.2 
1+0.4 

-0.5 
,+0.3 

-0.2 
-+0.4 

-0 
9 9+0.4 
'^•^-o.e 



0.73 
1-1^ 
4.6^ 
0.5"* 



1.22 
1.74 
1.51 
1.54; 
1.48 



+0.17 

0.17 
+0.28 

0.15 
+0.06 

0.06 
+0.13 

0.16 
+0.22 

0.12 



4.1 
1 

5.9 
2.3 
2.5 



+1.5 

0.8 
+0.3 
-0.4 
+0.5 

0.4 
+0.5 

0.3 
+0.4 

0.6 



0.50 
0.42 
2.04 
0.40 
0.37 



+0.05 

0.06 
+0.03 

0.03 
+0.05 

0.05 
+0.02 

0.03 
+0.03 

0.03 



0.5 

0.4 

2.3 

0.37 

0.37 



+0.5 

0.3 
+0.3 

0.2 
+0.5 

0.5 

+0.21 

0.17 

-fO.27 

0.19 



6.8/12 

27.0/28 
123.5/115 
37.4/26 
23.8/23 



Model: MCD + blackbody 



361+1302 
6097 
5644 
6361 

8190 



12 
5.3 
4.5 
7.0 
4.9 



+3. 
-2. 
+0.7 

0.6 
+0.4 

0.3 
+2.3 

1.1 
+0.7 

0.7 



^+0.10 

-0.08 

s+0.15 

-0.10 

-+0.014 

-0.013 

1 11+0-27 

3+0.15 
-0.14 



1.02: 

0.88J 
1.185'' 
l.lll 
1.59^ 



hi. 2 
-1.1 

1+0.08 
-0.42 
44+1-8 

,+0.4 
-0.4 



5.2: 

5.; 
5.323 

4.4"' 
2.1"! 



0.163 

0.21 

0.290 

0.20 

0.225 



+0.022 

0.019 
+0.03 
0.05 
+0.012 
0.027 
+0.06 
0.04 
+0.021 
0.018 



1+0.007 

-0.012 
-.+0.004 

-0.007 
^+0.009 

-0.016 
j+0.013 

-0.021 
n sq5+°°i2 



0.291] 
0.186J 
0.7641 
0.378J 



0.38 

0.21 

0.72 

0.4 

0.36 



+0.24 

0.18 
+0.23 

0.16 
+0.04 

0.11 
+0.5 
0.3 
+0.21 

0.17 



0.33l°:« 

+0.08 
-0.07 
1.8 
0.6 

+0.20 
0.11 



0.60 

i.o; 

0.48 



51.6/52 

36.5/48 

151.4/147 

17.5/23 

80.8/80 



10- 



Note. 

21 ^^- 



Col. (1): Observation ID; those in bold font are not in the low/hard state. Col. (2): Absorption column density 
Col. (3): Disk inner temperature in units of keV. Col. (4): Ri^ is the disk inner radius in units of 10 km and i 



10^^ ergs s ^. Col. (9): Best-fit x^ a-nd degrees of freedom. 



Col. (8): Bolometric luminosity of the blackbody component 



in units of 
is the disk 

— 2 —1 

cm s . 
in units of 



lar absorption on each spectrum from observations listed 
in Table [TJ The best-fit parameters are listed in Table [H 
The source luminosity was found at two different levels, 
a relatively low level of about (1 — 2) x 10"^^ ergs s~^ and 
a high level above 4 x 10'^^ ergs s~^. The long time scale 
lightcurve from all of the eighteen Chandra observations 
are shown in Figure [5] For the three short ACIS and the 
four HRC observations, which are inadequate for spectral 
fitting, their 0.3-8.0 keV luminosities are calculated using 
PIMMS from count rates, assuming a power-law spec- 
trum with the average absorption column density and 
power-law photon index. They are found to be consis- 
tent with the relatively low flux level. 

We then tried a multicolor accretion disk (MCD) 
model. A single MCD model with absorption did not 



improve the fits with respect to the power-law model. 
For five out of the ten spectra, an additional soft black- 
body component to the MCD was justified at a signifi- 
cance of 3 cr or more. Spectral parameters derived from 
the MCD model or the MCD plus blackbody are listed 
in Table [3l The bi-thermal model produced comparable 
or better fits than the power-law model. The disk lu- 
minosity versus disk inner temperature is plotted in the 
top panel of Figure [31 The soft blackbody component 
mainly lies in the low energy band, and the diskbb com- 
ponent contributes more than 90% of the unabsorbed 
fiux in the 2.0-8.0 keV band. We note that four of the 
five brightest observations require the soft component, 
but the brightest observation, 10027, does not. If we add 
a blackbody component into the 10027 spectrum and fix 
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Figure 2. Long term X-ray lightcurve of X37.8-I-54 from all avail- 
able Chandra observations. The diamond indicates the long ACIS 
observations and the luminosity is quoted from Table [2] The cross 
indicates short ACIS and HRC observations and the luminosity is 
estimated from detected count rates using PIMMS. 
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Figure 3. The disk bolometric luminosity (assuming face-on) 
versus the disk inner temperature of X37. 8-1-54 with values ob- 
tained from the best-fit model in Table [S] In the bottom panel, 
data points from the five observations in which the source flux is 
low and the spectrum is hard are removed. The dashed line indi- 
cates a best-fit power-law relation Ldisk oc T. ' 'to the data. 



the temperature at its mean value of 0.2 keV, the black- 
body luminosity is 0.5 x 10'^^ ergs s~^, similar to those 
from other observations, and it contributes only 2.3%, 
with a up limit of 5%, to the total observed flux. The 
MCD parameters are consistent with those from a single 
MCD model. In the other four spectra, the fraction of the 
soft blackbody is larger than 18%. The non-detection of 
the cool blackbody component in 10027 could be due to 
its relative weakness when the MCD component became 
much stronger than in others. 



Table 4 

Spectral parameters of X37. 8-1-54 from observation 10027 fitted 
with the KERRBB model subject to interstellar absorption. 



Parameters 



Values 



A/bh (Mq) 



m (IQi** g s-i) 

i (degree) 
xVdof 



5.0 ±0.3 

36 ±20 

16 2+^^-° 

0.9986 (> 0.80) 

52 ±16 
118.6/113 



Note. — Torque at the inner edge of the disk is zero, the 
hardening correction is 1.7, the distance to the source is 3.63 Mpc, 
and the self-irradiation and limb-darkening are turned on. 



Except for observations in which the source spectra 
are hard and the luminosities are the lowest, the disk 
luminosity seems to be correlated with the disk inner 
temperature. In the bottom panel of Figure [31 we re- 
moved points from observations when the source was 
in the low/hard state, specifically when Lx < 2 x 10"^^ 
ergs s~^ and F < 2.3, and then fitted the remaining ob- 
servations to a power-law function. There appears to be 
a robust correlation and the best-fit relation between the 
luminosity and temperature is Ldisk = T^^ ' with 1 a 
error. For the five observations (379, 2933, 8190, 10025, 
10026) in which the spectrum is hard and the luminos- 
ity is relatively low, the source may be in the hard state 
with a power-law spectrum. Due to low statistics, these 
spectra can also be fitted by thermal disk radiation, but 
it is not surprising that they do not vary in a Ldisk oc T^ 
pattern. 

The MCD model has been widely used as a standard 
indicator of the thermal dominant state. However, it 
does not include relativistic effects. We thus tried a more 
physical multicolor disk blackbody model with relativis- 
tic effects fully taken into account (KERRBB in XSPEC; 
iLi et al.ll2005[ ). to the spectrum from observation 10027 
when the source was brightest. In the model, we set zero 
torque at the inner boundary of the disk, a hardening 
correction of 1.7, a distance of 3.63 Mpc, and turned 
on self-irradiation and limb-darkening. The KERRBB 
model provides adequate fits to the spectrum with pa- 
rameters (absorption column density, black hole mass, 
mass accretion rate, spin, and inclination) listed in Ta- 
ble H 

Models consisting of a power-law with an additional 
soft blackbody or cool MCD component were also fitted 
to the spectra. Neither showed obvious improvement rel- 
ative to a single power-law model for most observations. 

For each observation, the short-term variability was 
examined using the Kolmogorov-Smirnov test. Variabil- 
ity in the 0.3-8 keV range was found only in observa- 
tions 5644 and 6361 with a significance of 9 and 5 a, 
respectively. In all other observations, the evidence for 
variability is at a confidence level below 3 a. We found 
that the source lies near the readout node boundary in 
observations 5644 and 6361, causing periodic oscillations 
when the source position on the detector plane moved 
across the node boundary due to dithering. The power 
spectrum shows a peak at 1000 s which is the dithering 
period along the direction perpendicular to the bound- 
ary. We therefore created a lightcurve from these two 
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Figure 4. Lightcurve of X37.8+54 in 0.3-8 keV created from in- 
tervals when the source is at least two pixels away from the readout 
node boundary. The horizontal bars indicate the intervals used for 
extracting energy spectra during the nondip and dip times, respec- 
tively. Simulations using MARX indicates that the variations seen 
in the figure are not from the instrument. 



observations with time steps of about 1000 s selected 
from intervals when the source was at least two pixels 
away from the readout boundary. Interestingly, there 
are dip-like variations on time scales of several ks. Sim- 
ulations using MARX 4.4.0 applying the aspect solution 
file and bad pixel map obtained from the observation 
indicates that dithering is not responsible for the vari- 
ations shown in Figure |4l We extracted energy spectra 
from intervals during the dip (58.171-65.171 ks) and the 
nondip (28.171-56.371 ks), which are marked in FigurelH 
We found that the dip spectrum is softer, with a power- 
law photon index of 3.5 ± 0.4 versus 2.65to'o9 from the 
nondip spectrum. The spectrum during dip is also softer 
than in the low state, which has a power-law index less 
than 2.3. Fitting the nondip spectrum with a blackbody 
plus MCD model modified by a shared absorption com- 
ponent leads to consistent results with from the entire 
spectrum. We then use two absorption components ap- 
plied onto each emission component individually. For the 
nondip spectrum, the absorption column density TVh on 
the MCD component is the same as on the blackbody 
component and is consistent with that derived using a 
single absorption component. For the dip spectrum, the 
fitting is made by fixing the temperature and A'h of the 
blackbody component to those of the nondip spectrum, 
due to limited photon number; other parameters, the 
normalization of the blackbody, the MCD and its ab- 
sorption are free in the fit. This leads to a best- fit A^h 
of 40.3^2^4^ X 10^^ cm^2 on the MCD component, which 
is significantly larger than the A'h from the nondip spec- 
trum, and suggests that the dip may be caused by extra 
absorption of the MCD component. The dip and nondip 
spectra with this blackbody plus MCD model are shown 
in Figure [5] for comparison. Observation 6361 was per- 
formed 62 ks after the end of 5644. The source flux in 
6361 seems to drop from the nondip level of 5644 to the 
relatively low level found in other observations. 

We searched in the literature and found no in- 
frared or radio counterpart for the X-ray source. 
The nearest radio source to X37. 8-1-54 is 37.54-1-53.2 
(|Rodriguez-rico et al.ll2004l) at an angular distance of 2" 
calculated using the X-ray position of R. A.=9^55'"46'^.61 
and decl.=-f69°4Q'41.1" (J20Q0.0) (.Kaaret et all [200l 



Figure 5. The dip (lower) and nondip (upper) spectra with their 
best fitted blackbody plus MCD model from observation 5644. 

iFeng fc Kaa"re^l2007bD . 

3. DISCUSSION 

We identify X37.8-h54 as the third ULX in the star- 
burst galaxy M82 based on its brightening in a recent 
Chandra observation to an observed fiux 5-10 times 
higher than previously observed. The peak isotropic lu- 
minosity in the 0.3-8 keV band corrected for interstellar 
absorption is (6.5 ± 0.3) x lO'^^ ergs s~^ inferred from 
fitting with a power-law model. 

The source was often found at relatively low luminos- 
ity, near (1 — 2) x 10'^^ ergs s~^. At the lowest lumi- 
nosities, < 1.5 X 10"^^ ergs s~^, the source appears to 
have hard spectra with the photon index less than 2.3. 
As the luminosity increases above this value, it shows 
dramatic spectral variation with the power-law photon 
index F varying from 1.9 to 3.0. This suggests a change 
in spectral state, and, indeed, most of the higher lumi- 
nosity spectra are better fitted with an MCD spectral 
model. 

The lower panel of Figure [3] shows that when observa- 
tions in the low/hard state are excluded, the remaining 
observations exhibit a robust correlation between Ldisk 
and Tin. The correlation coefficient between log L^isk and 
log Tin is 0.997 with a chance probability of 1.7 x 10~^. 
Fitting the data to a power-law relation, we obtained 
a power-law slope of 4.8 ± 0.9, consistent with the 4th- 
power relation predicted for the MCD model in the ther- 
mal dominant state. We note that the disk luminosity 
varies by about an order of magnitude. Thus, the data 
provide an adequate test of this relation. We conclude 
that X37. 8-1-54 exhibits two spectral states: a low/hard 
state and a thermal domin ant state. 

Following the recipe in iMakishima et al.l ()2000D and 
adopting canonical correction coefficients with k = 1.7 
and ^ = 0.412 , where k is the hardening correction factor 
and ^ is to correct for the radius of maximum temper- 
ature, the disk luminosity can be expressed as a func- 
tion of the black hole mass {m in solar mass) and the 
disk inner temperature (Tin in keV) as Ldisk = 7.2 x 
IQ'^^a^m^T^ ergs s~^, where a is related to the black 
hole spin with a — 1 for a non-spinning black hole and 
a — 1/6 for maximal spinning. Using a 4th power-law 
relation to fit the data points in the bottom panel of Fig- 
ured we obtain Ldisk = (3.7±0.3)I]y cosi lO^^ ergss"!. 
Thus, we have m = (7.2 ± 0.3)/(a-\/cosi), suggesting the 
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ULX contains a black hole of about 10 solar masses if 
it is a Schwartzchild black hole with a face-on disk, or a 
few tens solar masses if it is fast spinning or viewed at a 
high inclination. 

The best-fit compact object mass is 36 ± 20 M©, de- 
rived from fitting with the fully relativistic disk model, 
indicative of a stellar black hole slightly more massive 
than Galactic ones. We also tried different hardening 
correction factors from 1.5 to 1.9. The fitted mass varies 
by 20% with respect to that from k — 1.7. Interest- 
ingly, the specific angular momentum converges to its 
maximally allowed value, a, — 0.9986 with a lower limit 
of 0.80, suggesting the source contains a fast spinning 
black hole. The mass accretion rate in Eddington unit 
is about 1.0, which is larger than that found in the 
thermal d ominant state of Galactic BH Bs of no more 
than 0.6 (|McClintock fc RemillardI 120061 ). but is under 
the theoretical upper limit of about 10 for a stellar- 
mass black hole with a radiation domina ted accretion 
disk ()Begelmanll2"00^ . iHui fc KroliH (|2008[) investigated 
the spectra of six disk-dominated ULXs and found five 
sources had best-fit masses between 23 and 73 Mq. 
The three ULXs that have a derived mass > 25 Mq 
at a confidence > 99.9% also show rapid spin, which 
is quite similar to X37.8-f54. The best-fit black hole 
mass of X37. 8-1-54 is also similar with that of IC 10 X- 1 
(jPrestwich et al.ll2007l: [Silverman fc Filippenkdl2008[ ). It 
has been sugge sted that IC 10 X-1 h as a persistently high 
accretion rate (Barnard ct al. 2008); a rapidly spinning 
black hole is naturally expected since accretion would 
bring a fair amount of angular momentum to the com- 
pact object. 

We conclude that stellar Kerr black holes with a mass 
of several tens Mq could appear as ULXs with a ther- 
mal spectrum of a few keV in the 1-10 keV ba nd, as 
having been suggested bv lMakishima et al.l ( |2000( ). Such 
emission could also appear as a hard power-law spec- 
trum if the statistics is inadequate to detect the spec- 
tral curvature. These black holes are more massive than 
Galactic black holes, which are mostly around 10 Mq, 
and may be formed by the collapse of massive stars in a 
low-me tallicity environment and a ccount for a portion of 
ULXs ()Zampieri fc Roberts! [20091 ). 

Similar sources, Suzaku J1305-4931 in NGC 4945 and 
NGC 253 X-2, have been reported in the literature, 
which show a idisk oc T-^ correlation as thei r lumi- 
nosit ie s varied by a facto r of a bout 2-3 dlsobe et al.l 
120081: iKaiava fc PoutanenI |2009| ). From surveys of 
bright ULXs in nearby galaxies, many hard ULXs were 
found and interpreted as IMBHs in the low /hard state 
(jFeng fc Kaaretll2009t iKaaret fc Fe"^ [20091 ). For those 
sources, most of them present a constant spectral in- 
dex or hardness despite large variation of the luminosity. 
Therefore, those sources may have a different nature from 
X37.8-H54. 

It has been suggested that ULXs are a special 
case (face-on) of the Galactic microquasar SS 433, 
where super-Eddington accretio n drives a pow erful wind 
off the accretion dis k ( Fabrika fc Meschervakov 200!lt 
iBegelman et al.l[2006l ). The wind could b e optically thick 
and produce a b lackbody continuum ([King fc Pounds! 
120031 lKingjl2004D . This may be responsible for the soft 
excess observed in m any ULXs with a tempera ture of 
about 0.1-0.3 keV (Stobbart et al..,2006; Poutanen et al.l 



I2007t iFeng fc Kaa"retll2007a|) . When fitting spectra not 
from the low/hard state with an MCD component, all 
observations except 10027 need an additional blackbody 
component at a significance of at least 3 a. No obvi- 
ous correlation was found between the luminosity and 
temperature of the soft blackbody. However, the uncer- 
tainties are large and this result is not very constraining. 
The calculated surface area of the soft blackbody emis- 
sion is about 10^^ cm^ for all observations, corresponding 
to a spheric radius of 10^ cm, about 70i?g for a 10 so- 
lar mass black hole, where Rg is the gravitational radius 
equal to GM /(? . The absorption column density varies 
significantly between observations no matter what mod- 
els are used. Such large variation is likely internal to the 
binary, implying the outflow is inhomogeneous and time 
variable. We note that the soft blackbody appears at en- 
ergies where the absorption is important. Thus caution 
is warranted when interpreting it. 

Interesting timing behavior is seen in one of the ma- 
jor outbursts caught by the observations 5644 and 6361 
(Figure 131). The source shows soft dips on time scales of 
several ks. Soft dips could be caused by ab sorption of the 
inner accretion disk, like in 4U 1630—47 (jTomsick et al.l 
Il998| ). This interpretation is also consistent with our 
spectral analysis results in the penultimate paragraph 
of Section 2. We note that the soft blackbody compo- 
nent is assumed not to suffer extra absorption. A dense 
massive outflow from the disk could be the origin of the 
extra absorption and also the source of the soft black- 
body. Other mechanisms could also produce soft dips, 
e.g., the possible ejection of the inner part of accretion 
disk into relativistic jets in GRS 1915-1-105 ()Naik et al.l 
1200 It) . The time scale of ks of the dip in X37.8-h54 is 
longer than that in 4U 1630-47 and in GRS 1915+105 
by about an order of magnitude. More detailed study of 
the soft dip requires data with better quality. 

To summarize, X37. 8-1-54 provides us a good candidate 
of a stellar Kerr black hole being ultraluminous with a 
near Eddington accretion rate. This ULX could be repre- 
sentative of a class of ULXs with relatively low luminos- 
ity, i.e. a few 10"^^ ergs s~^. The luminosity of X37. 8-1-54 
is usually not high enough to be classified as a ULX, 
while in the outburst it truly enters the ultralumious 
regime. The X-ray spectrum at peak luminosity is domi- 
nated by a disk component, which is unlike many brighter 
ULXs at their maximum luminosity that they enter the 
'ultraluminous state' characterized by a st rong Comp- 
tonizing component ([Gladstone et al.l [20091 ). Therefore, 
this source may be in between the low luminosity tail of 
ULXs and high luminosity tail of Galactic BHBs. We 
note that the hard state of X37. 8-1-54 is more luminous 
than the canonical hard state found in Galactic BHBs. 
This is consistent with the correlation found in outbursts 
of BHBs that the peak luminosity of the hard state is 
scaled with the peak l uminosity of the thermal state in 
the following outburst (|Yu fc Yanll 20091 ). 
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